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Abstract 

A systematic theoretical study of the electronic and optical properties of Gai-^In^As 
self-assembled quantum-wires (QWR's) made of short-period superlattices (SPS) 

o : 

with strain-induced lateral ordering is presented. The theory is based on the effec- 
tive bond-orbital model (EBOM) combined with a valence-force field (VFF) model. 
Valence-band anisotropy, band mixing, and effects due to local strain distribution at 

m ' 

I/-) . the atomistic level are all taken into account. Several structure models with varying 

o , 

degrees of alloy mixing for lateral modulation are considered. A valence force field 

o 

model is used to find the equilibrium atomic positions in the QWR structure by min- 

!i2 i 

imizing the lattice energy. The strain tensor at each atomic (In or Ga) site is then 
obtained and included in the calculation of electronic states and optical properties. 

o 
o 



X 



It is found that different local arrangement of atoms leads to very different strain 
distribution, which in turn alters the optical properties. In particular, we found that 
in model structures with thick capping layer the electron and hole are confined in the 
Ga-rich region and the optical anisotropy can be reversed due to the variation of lat- 
eral alloying mixing, while for model structures with thin capping layer the electron 
and hole are confined in the In-rich region, and the optical anisotropy is much less 
sensitive to the lateral alloy mixing. 
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I. INTRODUCTION 



Optical properties of III-V semiconductor nanostructures have attracted a great deal of inter- 
est recently for their applications in optical communications that involve switching, amplification, 
and signal processings. Ga x In\_ x As and Ga x Irii_ x P are among the most important ternary III-V 
compound semiconductors. The band gap of Ga x Ini_ x As covers both the 1.3 and 1.55 /an wave- 
lengths, which are preferred in long distance fiber communications. However, long- wavelength pho- 
tonic devices based on lattice-matched Gao.47ln 0j53 As/InP heterostructures suffer from strong Auger 
recombination and inter- valence band absorption processes [1-3]. Recently, to improve the perfor- 
mance of long- wavelength semiconductor lasers, long wavelength(~ 1.55/zm) Ga^In^s-As quantum- 
wire (QWR) lasers have been grown by a single step molecular beam epitaxy technique [4-6]. 

An important optical property is the change in the emission of light (in energy, polarization, 
and intensity) that results from phase-space filling of carriers in one- and two-dimensionally (2D) 
confined systems. As the dimensionality of the quantum confinement increases from ID to 2D, 
the sharpening of the density of states gives rise to a lower excitation threshold, thereby yielding 
potentially enhanced optical effects. It is found that the QWR laser structure is a promising candidate 
for optical communication device because of the many predicted benefits, such as higher gain, reduced 
temperature sensitivity, higher modulation bandwidths, and narrower spectral linewidths[5]. 

The fabrication of quantum wires via the strain-induced lateral-layer ordering (SILO) process 
starts with the growth of short-period superlattices (SPS) [e.g. (GaAs)2/(InAs)2.25] along the [001] 
direction. The excess fractional InAs layer leads to stripe-like islands during the MBE growth. [4] 
The presence of stripes combined with the strain effect lead to a natural phase separation as ad- 
ditional layers of GaAs or InAs are deposited. A self-assembled quantum wire heterostructure can 
then be created by sandwiching the composition modulated layer between barrier materials such as 
Alo.24Gao.24Ino.52 As (quarternary), Alo. 48^0.52 As(ternary), or InP (binary) [4-6] . It was found that 
different barrier materials can lead to different degree of lateral composition modulation, and the 
consequent optical properties are different. [6] Besides the self-assembled lateral ordering, it is be- 
lieved that the strain also plays a key role [5,6] in the temperature stability and optical anisotropy for 
the QWR laser structure. Much work has been undertaken which predicts that by using strained- 
layer superlattice to form the active region of a quantum-wire laser, the threshold current can be 
decreased by one order of magnitude, and the optical loss due to intervalence-band absorption and 
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Auger recombination will also be greatly reduced [4- 6]. Also, the temperature sensitivity is reduced 
by an order of magnitude compared with strain-free structures. Temperature sensitivity of the las- 
ing wavelength for typical GalnAs/InP lasers is around 5A/°C. By using a distributed-feedback 
structure, the temperature dependence of the lasing wavelength can be reduced to lA/°C[5]. With 
the self-assembled GalnAs QWR, the temperature sensitivity is smaller than 0.lA/°C[5,6]. Thus, it 
represents an important improvement on the current technology in fabricating long wavelength lasers 
for fiber communication. 

A few theoretical studies of Ga^Ini-^As self-assembled QWRs based on a simplified uniform strain 
model have been reported. [8-10] In these calculations the SPS region is modeled by a Ga x Ini_ x As 
alloy with a lateral modulation of the composition x. Although these calculations can explain the 
QWR band gap and optical anisotropy qualitatively, it does not take into account the detailed SPS 
structure and the microscopic strain distribution. The understanding of these effects is important 
if one wish to have a full design capability of the self-assembled QWR optoelectronic devices. On 
the other hand, microscopic strain distributions in several SPS structures have been calculated, and 
the electronic properties of these SPS structures have been studied via an empirical pseudopotential 
method[ll]. However, the effects of sandwiching the SPS structures between barrier materials so 
as to form QWRs have not been explored. It is also worth pointing out that the effective masses 
obtained in the pseudopotential method for the constituent bulk materials are 0.032m and 0.092m 
for InAs and GaAs, respectively. [12] These values are 30% larger than the actual values (0.024m and 
0.067m ); thus, the energy levels of quantum confined states obtained by this method are subject to 
similar uncertainty. 

In this paper we present a systematic study of the effects of microscopic strain distribution on 
the electronic and optical properties of the Gai-^In^As self-assembled QWRs via the combination 
of the effective bond-orbital model (EBOM) for electronic states and the valence-force-field (VFF) 
model for microscopic strain distribution. Both clamped and undamped SPS structures with different 
degrees of lateral alloy mixing are considered. The clamped structure has a SPS layer sandwiched 
between the substrate and a thick capping layer so that the interface between the top monolayer 
of the SPS and the barrier region is atomically flat. The undamped structures correspond to a 
SPS layer sandwiched between the substrate and a thin capping layer, which allows more flexible 
relaxation of interface atoms so the strain energy can be relieved. This leads to a wavy interface 
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structure. Most self-assembled QWR samples reported to date are closer to undamped structures 
with a wavy interface, although it should be possible to produce self-assembled QWRs closer to the 
clamped structures by using a thick capping layer. 

Our theoretical studies show that there are profound differences in the electronic and optical 
properties between the clamped and undamped self-assembled QWR structures. In particular, we 
found that in clamped structures the electron and hole are confined in the Ga-rich region and the 
polarization of Photoluminescence (PL) can go from predominantly along [110] for SPS with abrupt 
change in In/Ga composition to predominantly along [110] for SPS with smooth change in In/Ga 
composition. On the other hand, for undamped structures the electron and hole are confined in the In- 
rich region, and the optical polarization is always predominately along [110] with a weak dependence 
on the lateral alloy mixing. This also implies that by changing the degree of strain relaxation at 
the interface between the SPS and the capping layer, one can tailor the optical properties of self- 
assembled QWRs. Thus, through this theoretical study, we can gain a better understanding of the 
strain engineering of self-assembled QWR structures which may find applications in fiber-optical 
communication. 

II. MODEL STRUCTURES 

For both clamped and undamped SPS structures, we consider different degrees of lateral alloy 
mixing and examine the effects of the microscopic strain distribution on the electronic and optical 
properties. Two example QWR model structures considered in the present paper (prior to alloy 
mixing) are depicted in Figure 1. The supercell of the first model structure consists of 8 pairs 
of (001) (GaAs) 2 /(InAs) 2 SPS with a total thickness of ps 9AA (quantum well region) followed by a 
Alo.24Gao.24Ino.52As layer (barrier region) with thickness ~ 60 A (20 diatomic layers). The supercell of 
the second model structure consists of 8 pairs of (001) (GaAs) 2 /(InAs)2. 2 5 SPS with a total thickness 
of ps 100 A (quantum well region) followed by a Alo.24Gao.24Ino.52As layer ( barrier region). We 
assume an arrangement of alternating stripe like islands due to strain induced lateral ordering. In 
the diagram, no lateral alloy mixing is shown. In our calculations of strain distribution and electronic 
structures, different degrees of lateral alloy mixing will be considered. Experimentally, self-assembled 
InGaAs QWRs were usually grown with the (2/2.25) SPS structure. [12] However, the migration 
of excess In atoms during MBE growth could lead a structure somewhere between (2/2) SPS and 
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(2/2.25) SPS structures. 

In both model structures, we can divide the self-assembled QWR into two regions with the left 
half being Ga rich and the right half In rich. During growth, varying degrees of lateral alloy mixing 
of these islands with the surround atoms is likely to occur. In the atomic layers with In-rich alloy 
filled the right half of the unit cell (layer 2 in structure 1 and layers 7 and 9 in structure 2), the In 
composition xj n is assumed to vary in the [110] direction (or y' direction) according to the relation 

x m [l - sm(ny'/2b)]/2 ioiy'<b 

for b < y' < L/2 - b 

•<■/» { x m {l + sin[7r(y' - L/2)/26]}/2 for L/2 - b < y' < L/2 + b (1) 

x m for L/2 + b<y'<L-b 

x m {l - sin[7r(y' - L)/26]}/2 for y> > L - b, 

where x m is the maximum In composition in the layer, 2b denotes the width of lateral composition 
grading, and L is the period of the lateral modulation in the [110] direction. In layer 4 of structure 1, 
which contains Ga-rich alloy in the left half of the unit cell, the Ga composition (xgo) as a function 
of y' is given by a similar equation with the sign of the sine function reversed as compared to Eq. 
(1). In structure 2, there are a few atomic layers that contain 0.25 monolayer of In (layers 3 and 13) 
or Ga( layers 5 and 11). The lateral alloy modulation in layers 3 and 13 is described by 

for < y' < 5L/8 - b 

x m {l + sm[n(y' - 5L/8)/26]}/2 for 5L/8 - b < y' < 5L/8 + b, 

x m for 5L/8 + b < y' < 7L/8 - b 

x m {l - sm[n(y' - 7L/8)/2b]}/2 for 7L/8 - b < y' < 7L/8 + b 

for 7L/8 + b<y' < L. 



(2) 



Similar equation for xc a in layers 5 and 11 can be deduced from the above. By varying the parameters 
x rn and b, we can get different degrees of lateral alloy mixing. Typically x m is between 0.6 and 1, 
and b is between zero and 15a[no] ~ 62A. 

A valence force field (VFF) model[13-15] is used to find the equilibrium atomic positions in the 
self-assembled QWR structure by minimizing the lattice energy. The strain tensor at each atomic 
(In or Ga) site is then obtained by calculating the local distortion of chemical bonds. We found 
that different local arrangement of atoms can lead to very different strain distribution. In particular, 



the shear strain in the clamped structure can change substantially when the In/Ga composition 
modulation is changed. Consequently, the optical anisotropy can be reversed due to the change in 
the strength of the shear strain caused by the intermixing of In and Ga atoms. 

III. THEORETICAL APPROACH 

The method used in this paper for calculating the strained QWR band structure is based on 
the effective bond-orbital model (EBOM). A detailed description of this method has been published 
elsewhere [16, 17]. EBOM is a tight binding-like model with minimum set of localized basis (the 
bonding or anti-bonding orbitals). The interaction and optical parameters are obtained by a cor- 
respondence with the k • p theory, which can be cast into analytic forms. Thus, the model can be 
viewed as a spatially discretized version of the k • p method, while retaining the virtues of LCAO 
(linear combination of atomic orbitals) method. The k • p model is the most popular one for treating 
electronic structures of semiconductor quantum wells or superlattices. However, when applied to 
complex structures such as self-assembled quantum wires[4-9] or quantum dots[13,18,19], the method 
becomes very cumbersome if one wish to implement the correct boundary conditions that take into 
account the differences in k • p band parameters for different materials involved. EBOM is free of 
this problem, since different material parameters are used at different atomic sites in a natural way. 
For simple structures, when both EBOM and k • p model are equally applicable, the results obtained 
are essentially identical. [16] 

The optical matrix elements for the QWR states are computed in terms of elementary optical ma- 
trix elements between the valence-band bond orbitals and the conduction-band orbitals. The present 
calculation includes the coupling of the four spin-3/2 valence bands and the two spin- 1/2 conduction 
bands closest to the band edges. Thus, it is equivalent to a 6-band k ■ p model. For our systems 
studied here, the band-edge properties are relatively unaffected by the split-off band due to the large 
spin-orbit splitting as discussed in our previous paper [9]. Hence, the split-off bands are ignored here. 
The bond-orbitals for the GaAs and InAs needed in the expansion of the superlattice states contain 
the following: four valence-band bond orbitals per bulk unit cell, which are p-like orbitals coupled 
with the spin to form orbitals with total angular momentum J=3/2 plus two conduction-band bond 
orbitals with J=l/2. They are written as 



\R,uj M >= ^C(a,(T, J,M)\R,a > tp a , (3) 

a,cr 

where J = 1/2 and 3/2 for the conduction and valence bands, respectively, and M = — J, • • • J, ^V, 
designates the electron spinor(<r=l/2,-l/2), and |R, a > denote an a-like(o; = s, x, y, z) bond orbital, 
located at unit cell R. C(a,a, J, M) are the coupling coefficients obtainable by group theory. All 
these bond orbitals are assumed to be sufficiently localized so that the interaction between orbitals 
separated farther than the nearest-neighbor distance can be ignored. 

The effect of strain is included by adding a strain Hamiltonian H st to the EBOM Hamiltonian[17]. 
The matrix elements of H st in the bond-orbital basis can be obtained by the deformation-potential 
theory of Bir and Pikus[20]. We use the valence-force field (VFF) model of Keating[14] and Martin[15] 
to calculate the microscopic strain distribution. This model has been shown to be successful in fitting 
and predicting the elastic constants of elastic continuum theory, calculating strain distribution in a 
quantum well, and determining the atomic structure of III-V alloys. It was also used in the calculation 
of the strain distribution in self-assembled quantum dots. [13,19] The VFF model is a microscopic 
theory which includes bond stretching and bond bending, and avoids the potential failure of elastic 
continuum theory in the atomically thin limit. The total energy of the lattice is taken as 

13 13^^ 
y = ^ X] ~^ a ij(dij — d 0ij ) I d 0ij + - ^2 ^ -^(3ijk{dij-dik + d 0ti jd 0tik /3) / d 0t ijd 0tik (4) 

ij i jy^k 

where % runs over all the atomic sites, j, k run over the nearest-neighbor sites of i, dij is the vector 
joining the sites i and j, d^ is the length of the bond, do,ij is the corresponding equilibrium length 
in the binary constituents, and and are the bond stretching and bond bending constants, 
respectively, a and (3 are from Martin's calculations. [11] For the bond-bending parameter (3 of In- 
As-Ga, = y/p[ij]P[ik] following Ref. [15]. 

To find the strain tensor in the InAs/GaAs self- assembled QWR, we start from ideal atomic 
positions and minimize the system energy using the Hamiltonian given above. Minimization of the 
total energy requires one to solve a set of coupled equations with 3N variables, where N is the total 
number of atoms. Direct solution of these equations is impractical in our case, since the system 
contains more than 6,000 atoms. We use an approach taken by several authors[10,13,19] which has 
been shown to be quite efficient. In the beginning of the simulation all the atoms are placed on the InP 
lattice, we allow atoms to deviate from this starting positions and use periodic boundary conditions 
in the plane perpendicular to the growth direction, while keeping atoms in the planes outside the 
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SPS region at their ideal atomic positions for a InP lattice (since the self-assembled QWR is grown 
epitaxially on the InP substrate). In each iteration, only one atomic position is displaced and other 
atom positions are held fixed. The direction of the displacement of atom % is determined according 
to the force = —dV/dxi acting on it. All atoms are displaced in sequence, the whole sequence 
is repeated until the forces acting on all atoms become zero at which point the system energy is a 
minimum. Once the positions of all the atoms are known, the strain distribution is obtained through 
the strain tensor calculated according to the method described in Ref. [21]. 
Let R° be the position matrix without strain: 



R° = 



( pO pO pO \ 

rl 12x -""232: rl 34x 

DO pO pO 

rl 12j/ ""-232/ ""-342/ 

pO pO pO 

Y - fl 12z -"-232 rl 342 J 



R° — R°;?, j = 1,4 and R° (i = 1,4) denote positions of four As atoms surrounding a 



where R°- 

Ga or In atom. Here we choose R° 



12 



;i,-l,0)a/2,R: 
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1, 0, -l)a/2, and K° u = (1, 1, 0)a/2. a 



is the lattice constant of GaAs or InAs, depending on the site. Let R be the corresponding position 
matrix with strain. It was shown in Ref. [21] that the strain tensor is given by 



R * Rq 1 - 1. 



(5) 



After getting the strain tensor, the strain Hamiltonian is given by Bir and Pikus[20] 



( -AV H + D 1 V3d 



St 



H 



where — (e^ + €ji)/2, and 



V 



xy 



xy 

AV H + D 2 



y/3de 



\ 



yz 



^yz 



-AV H + D 3 



(6) 



AV H — (a 1 +a 2 )(£ xx +e y y+e zz ), D x — b(2e xx — e yy — e zz ), D 2 — b(2e yy — e xx — e 22 ), D 3 — b(2e zz — e yy — e xx ), 
The strain potential on the s states is given by 

AV C = ci(e xx + e yy + e zz ), 



The strain Hamiltonian in the bond-orbital basis \ JM > can be easily found by using the coupling 
constants[13], i.e, 



< JM\H st \J'M' >= C ( a , o\ J > M)*C(a', a; J', M')H s ^ a , (7) 

The elastic constants C± 2 and Cn for GaAs, InAs and AlAs can be found in Ref. [22,23]. The 
deformation potentials a±, a 2 , b, c±, d can be found in Ref. [23-26]. The linear interpolation and 
virtual crystal approximations used to obtain the corresponding parameters for the barrier material 
(Alo.24Gao.24Ino.52 As). 

The above strain Hamiltonian is derived locally for the each cation atom in the self-assembled 
QWR considered. To calculate the electronic states of the self-assembled QWR for model structure 
l[see Fig. 1(a)], we first construct a zeroth-order Hamiltonian for a superlattice structure which 
contains in each period 8 pairs of (001) (GaAs)2/(InAs)2 SPS layers (with a total thickness around 
100 A) and 20 diatomic layers of Al .24Gao. 24I110.52 As (with thickness around 60 A). So, the superlattice 
unit cell for the zeroth-order model contains 52 diatomic layers. The appropriate strain Hamiltonian 
for the the (GaAs) 2 /(InAs) 2 SPS on InP is also included. For model structure 2 [see Fig. 1(b)], the 
zero-order superlattice consistes of two additional monolayers of InAs inserted into the 8 pairs of 
(2/2) SPS: one between the 2nd and 3rd pair of (2/2) SPS, the other between the 5th and 6th pair 
of (2/2) SPS. 

The eigen-states for the zero-th order Hamiltonian for different values of k 2 (separated by the SL 
reciprocal lattice vectors in the [110] direction) are then used as the basis for calculating the self- 
assembled QWR electronic states. The difference in the Hamiltonian (including strain effects) caused 
by the intermixing of Ga and In atoms at the interfaces is then added to the zeroth-order Hamiltonian, 
and the electronic states of the full Hamiltonian is solved by diagonalizing the Hamiltonian matrix 
defined within a truncated set of eigen-states of the zeroth-order Hamiltonian. A total of ~ 500 eigen- 
states of the zero-th order Hamiltonian (with 21 different k nQ points) were used in the expansion. 
The subbands closest to the band edge are converged to within 0.1 meV. 

IV. RESULTS AND DISCUSSIONS 

A. Strain distributions 

In this section we discuss strain distributions in our model structures as described in section I 
with varying degree of lateral alloy mixing. Both clamped (with thick capping layer) and undamped 
(with thin capping layer) situations are considered. To model the alloy structure with composition 



modulation, we use a super-cell which contains 72 atoms in the [110] {%)') direction, 36 atoms in 
the [110] (x r ) direction and 64 or 68 atomic planes along the [001] (z) direction [8 pairs of (2/2) or 
(2/2.25) SPS] plus a GaAs capping layer. In the atomic planes which consist of alloy structure, we 
first determine the In composition at a given y' according to Eq. (1) and then use a random number 
generator to determine the atomic species along the x' direction. The bottom layer of atoms are 
bonded to the InP substrate with the substrate atoms fixed at their ideal atomic positions. The 
calculated strain distributions are then averaged over the x' coordinate. For the clamped case, the 
GaAs capping layer is assumed to be lattice matched to InP with a flat surface. For the undamped 
case, the capping layer is allowed to relax freely, thus giving rise to a wavy surface structure. 

The diagonal strains in the four atomic layers that constitute the (2/2) SPS in structure 1 for the 
clamped and undamped cases are shown in Figs. 2 and 3, respectively. The lateral alloy modulation 
considered is described by Eq. (1) with x m = 0.8 and b = 7a[no]. For best illustration, we show 
diagonal strains in a rotated frame, in which x' is [110], y' is [110], and z' is [001]. The layer number 
in the figure labels the atomic layers in Fig. 1(a), starting from the bottom layer as layer 1. There are 
two main features of the microscopic strain distribution. First, in the ideal situation (without atomic 
relaxation), one would predict e y i y i to be the same as e x > x i due to symmetry. However, with atomic 
relaxation, all Ga(In) atoms tend to shift in a direction so as to reduce the strain in the Ga(In)-rich 
region. Thus, the magnitude of e y > y > on Ga(In) sites (dashed lines) is lower(higher) than e x r x > in 
Ga(In)-rich region. Second, the z component strain (dash-dotted lines) tend to compensate the other 
two components such that the volume of each bulk unit cell is closer to that for the unstrained bulk. 
Thus, we see that e z v has an opposite sign compared to e x > x i or e y / y i at all atomic sites. 

Comparing Fig. 2 with Fig. 3, we see that the main difference between the clamped and un- 
damped case is that in the atomic layers with lateral composition modulation (layers 2 and 4), the 
hydrostatic strain (sum of all three diagonal strain components) is much smaller in the undamped 
structure than in the clamped structure. This would lead to a major difference in the band-edge 
profile for the two structures to be discussed below. 

The difference in e x / x > and e y > y > leads to a nonzero shear strain e xy = (e y > y > — e x > x ')/2 in the original 
coordinates. For the clamped case, the shear strain is found to be particularly strong near the 
boundary where the alloy composition begins to change, and it is sensitive to the degree of lateral 
alloy mixing. For x m = 0.8 and b = 7a[no], the maximum value of e xy is around 0.4% For abrupt 
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composition modulation (x m = 1 and 6 = 0) the maximum e xy value increases five-fold to around 
2%. The other shear strain components (e xz and e yz ) are found to have similar magnitude. For the 
undamped case, the shear strain is strong even in regions away from the boundary, and the maximum 
shear strain is larger than its counterpart in the clamped structure by about 30%. 

The diagonal strains of structure 2 for the undamped case with x m — 1 and b = 7a[ 110 ] are shown 
in Fig. 4. The layer number in the figure labels the atomic layers in Fig. 1(b), starting from the 
bottom layer as layer 1. Only four representative atomic layers (3,4,5, and 7) are shown. We note 
that the average magnitude of the hydrostatic strain in the In-rich region in model structure 2 is 
comparable to that for undamped structure 1. 

B. Electronic structures 

In order to understand the aspect of lateral quantum confinement due to composition modulation 
and strain, we examine the band-edge energies of a strained quantum well structure whose well 
material is the same structure as appeared in the SILO QWR with a fixed value of y' and the barrier 
is 60Athick Alo.24Gao.24Ino.52As layer. For structure 1 depicted in Fig. 1(a), the well material consists 
of 8 pairs of (GaAs) 2 /(InAs) 2 SPS. The conduction band minimum and valence band maximum of 
the above quantum well with different degrees of lateral alloy mixing as functions of y' (dashed 
lines) are shown in Fig. 5 for both clamped and undamped cases. The strain Hamiltonian used 
here is the same as the one used in the self-assembled QWR at the corresponding y'. All material 
parameters are chosen the same as in Ref. [9] for temperature at 77K, except that the deformation 
potential c\ of GaAs is slightly modified from — 6.8eV to — 7.1eV so that ai + a 2 + c\ = — 9.8eV 
is in agreement with the experimental measurements [26] and the valence-band offset between GaAS 
and InAs used here is 0.26 eV according to the model solid theory[27]. To correct the band gap of 
SPS and alloy due to the bowing effect, we add a correction term — 0.4a;(l — x)eV to the diagonal 
element for the s-like bond orbital in the Hamiltonian, where x is the effective alloy composition 
of the SPS or alloy at a given y'. With the correction of bowing effect, our model gives band 
gaps for the (GaAs) 2 (InAs) 2 SPS, Gai-^In^As alloy, and the superlattice made of (GaAs) 2 (InAs) 2 
SPS grown on InP substrate all in very good agreement with available experimental observations. 
The photoluminescence (PL) measurements indicates that the (GaAs)2(InAs)2 SPS grown on InP 
substrate has a gap around 0.75 eV.[28] Here we obtain a band gap of 0.74 eV for the (GaAs) 2 (InAs) 2 
SPS (with a strain distribution obtained again via the VFF model) and 0.79 eV for 8 pairs (or 100 
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A) of (GaAs) 2 (InAs) 2 SPS sandwiched between Alo.24Gao.24Ino.52As confining barriers. This is also 
consistent with the PL measurements on the (GaAs) 2 (InAs) 2 /InP multiple quantum wells. [28] 

The band-edge profiles shown in this figure suggest that for the clamped case with alloy mixing, 
both electrons and holes are confined in the Ga-rich region with an effective barrier height around 
0.13eV for the electron and O.leV for the hole. Both offsets are large enough to give rise to strong 
lateral confinement for electrons and holes in the Ga-rich region. Without alloy mixing, the electron 
is not well confined, since the average energy between the Ga-rich and In-rich regions are nearly the 
same. For the undamped case, both electrons and holes are confined in the In-rich region with an 
effective barrier height around 0.2 — 0.3eV for the electron and 0.1 — 0.2eV for the hole. We note 
that the band gap in the In-rich region is rather insensitive to the alloy mixing, while the effective 
barrier height can change substantially due to different alloy mixing. 

Fig. 6 shows the band-edge profiles for 8 pairs of (2/2.25) SPS [as illustrated in Fig. 1(b)] 
sandwiched between 6OAAlo.24Gao.24Ino.52As barriers with and without alloy mixing for both clamped 
and undamped cases. For the clamped case, the electron is confined in the In-rich region, while the 
hole is confined in the Ga-rich region. Thus, we have a spatially indirect QWR, which will have very 
weak inter-band optical transition. For the undamped case, the band-edge profile is similar to the 
undamped case of structure 1 with both the electron and hole confined in the In-rich region. The 
band gaps for all cases of structure 2 are consistently smaller than their counterparts in structure 1. 

Fig. 7 shows the near zone-center valence subband structures and squared optical matrix elements 
(P 2 ) for the QWR model structure 1 [as depicted in Fig. 1(a)] with thick capping layer (clamped case) 
and without alloy mixing. Here P 2 is defined as J2s,s> I < tyv\£ • plV'c > | 2 > where i\) v (ip c ) denotes a 
valance (conduction) subband state. The symbol J2 S s ' means a sum over two near degenerate pair of 
subbands in the initial and final states, e denotes the polarization of light. Here we consider only the 
x' (along the QWR axis) and y' (perpendicular to the QWR axis) components. Only the dispersion 
along the k\ ([110]) direction is shown, since the dispersion along the k 2 direction for the confined 
levels is is rather small due to strong lateral confinement. All subbands are two- fold degenerate at 
the zone center due to the Kramer's degeneracy and they split at finite wave vectors as a result of 
lack of inversion symmetry in the system. The first three pairs of subbands are labeled VI, V2, and 
V3. They have unusually large energy separations compared with other valence subbands. This is 
because the first three pairs of subbands represent QWR confined states, whereas the other valence 
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subbands (with energies below -110 meV) are unconfined by the lateral composition modulation. 

To examine the effects caused by the shear strain, we have also calculated the band structures 
with the shear strain set to zero. We found that without the shear strain the first pairs of valence 
subbands (VI) have much larger dispersion along the k\ direction with an effective masses along [110] 
about a factor five smaller than that with the shear strain. This indicates that the VI subband states 
for the case without shear strain are derived mainly from bond orbitals with ^'-character, which leads 
to stronger overlap between two bond orbitals along the x' direction, hence larger dispersion along 
k\. When the shear strain is present, the character of bond orbitals in the VI subband states change 
from mainly x'-like to predominantly y'-like; thus, the dispersion along k\ becomes much weaker. 
This explains the very flat VI subbands as shown in Fig. 6. The switching of orbital character in the 
VI subbands can be understood as follows. When the shear strain is absent, we have e x ' x > = e y > y >, and 
the confinement effect in the y' direction pushes down the states with y' character (which has smaller 
effective mass in the y' direction), thus leaving the top valence subband (VI) to have predominantly 
x' character. On the other hand, with the presence of shear strain as shown in structure 1, we have 
€yi y i < e x i x i in Ga-rich region, and the bi-axial term of the strain potential forces the y'-like states 
to move above the a/- like states, overcoming the confinement effect. As a result, the VI subbands 
become y'-like. The conduction subbands are approximately parabolic as usual with a zone-center 
subband minimum equal to 688 meV This gives an energy gap of 763 meV for QWR model structure 
1. As seen in this figure, the Cl-Vl transition for ^/-polarization (dashed line) is more than twice that 
for the a:'-polarization (solid line) with -P[ii ]/-P[ii ] ~ 0-36. This is consistent with the discussion for 
valence subband structures, where we concluded that the bond orbitals involved in the VI subbands 
are predominantly j/'-like. 

Next, we study the case with alloy mixing, in which a gradual lateral modulation in In composition 
substantially reduces the shear strain. Using a lateral alloy mixing described in Eq. (1) with x m = 0.8 
and b = 7a[n ], we found significant change in band structures and optical properties compared to 
the case without alloy mixing. With alloy mixing we found that the first pairs of valence subbands 
(VI) have much larger dispersion than their counterparts in Fig. 6, and the atomic character in the 
VI states change from predominantly y'-like (without alloy mixing) to rr'-like (with alloy mixing). 
The change of atomic character in the VI states is a consequence of the competition between the 
QWR confinement (which suppresses the y' character) and the e xy shear strain (which suppresses the 
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x' character). As a result, the ratio P[ii ]/P[ii ] changes from 0.36 (without alloy mixing) to 3.5 (with 
alloy mixing), indicating a reversal of optical anisotropy. The above studies lead to the conclusion 
that for self-assembled QWRs with thick capping layer (clamped case), the optical anisotropy is 
sensitive to the degree of lateral alloy mixing. 

Next, we study the properties of self-assembled QWRs with thin capping layer (undamped case). 
We found that for the undamped case, the electron and hole are always confined in the In-rich 
region, regardless of the degree of alloy mixing and whether the QWR structure consists of (2/2) 
SPS or (2/2.25) SPS. Furthermore, since the strain in In-rich region has opposite sign compared 
to that in the Ga-rich region, the effect of shear strain on the VB states confined in In-rich region 
is also reversed. Namely, it pushes the energy of the y'-like VB states down relative to the x'-like 
VB states. Thus, both the QWR confinement and the shear strain effects are in favor of having a 
predominant x' character in the VI subband. Therefore, the inter-band optical transition always has 
a ratio P[ii ] /-P[no] larger than 1. 

To illustrate this, we show in Fig. 8 the near zone-center valence subband structures and squared 
optical matrix elements (P 2 ) for the QWR model structure 1 [as depicted in Fig. 1(a)] with thin 
capping layer (undamped case) and with alloy mixing described by x m = 0.8 and b = 7a[n ]. 
Comparing the band structure in Fig. 8 with that in Fig. 7, we see that the VI subband now has 
much larger dispersion (i.e. much smaller effective mass) than its counterpart in Fig. 6. There are 
two reasons for this. First, the hole in the undamped case is now confined in the In-rich region rather 
than in the Ga-rich region as in the clamped case, thus having much smaller effective mass. Second, 
both the QWR confinement effect and the shear strain effect lead to a predominantly x' character in 
the VI states, which tend to yield a smaller effective mass in the k\ direction. The optical properties 
shown in Fig. 8 are also distinctly different from that shown in Fig. 7. The inter-band optical 
transition (Cl-Vl) now has a much stronger polarization along [110] than [110]. The quick drop of 
the the y' polarization in Cl-Vl transition at k\ 0.035^ is caused by the band mixing of VI and 
V2 states as can be seen in the VB band structures in this figure. 

We have also calculated the near zone-center valence subband structures and squared optical 
matrix elements (P 2 ) for the QWR model structure 2 with thin capping layer (not shown). We 
found that the band structures and optical matrix elements are similar to Fig. 8, except that the 
energy separation between VI and V2 subbands are smaller (by ~ 5mV) and the polarization ratio 
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( P [ 2 iio]/ P [no]) is lar S er - 

Finally, we list in Table I the band gaps and the polarization ratio (-P[iioi/-P[iio]) f° r var i°us QWR 
structures with different degrees of lateral alloy mixing (as indicated by the parameters x m and b) . 
To calculate the polarization ratio in the photoluminescence (PL) spectra, we integrate the squared 
optical matrix element over the range of k\ corresponding to the spread of exciton envelope function 
in the k\ space. The exciton envelope function is obtained by solving the ID Schrodinger equation 
for the exciton in the effective- mass approximation similar to what we did in Ref. 9. For undamped 
(2/2) SPS structure, the band gap varies between 0.77 and 0.8 eV and the polarization ratio changes 
from 3.1 to 1.5 as the degree of alloy mixing is increased. For undamped (2/2.25) SPS structure, 
both the band gap (around 0.74 eV) and the polarization ratio (around 2) are insensitive to the 
degree of alloy mixing. 

Experimentally, the PL peak for self-assembled InGaAs QWRs with similar dimensions as con- 
sidered in our model structures are around 735 meV and the polarization ratio for different samples 
are found to be between 2 and 4. [4-6] Taking into account an exciton binding energy between 10 
and 20 meV, we find that our theoretical predictions for both undamped and clamped structures 
with lateral alloy mixing are consistent with the experimental finding. To sort out which one is 
closer to the experimental structure would require more detailed comparison between the theory and 
experiment for optical spectra involving the excited hole states such as the PL excitation spectra. 

V. CONCLUSION 

We have calculated the band structures and optical matrix elements for the self-assembled GalnAs 
QWR grown by the SILO method. The actual SPS structure and the microscopic strain distribution 
has been taken into account. The effects of microscopic strain distribution on the valence subband 
structures and optical matrix elements are studied for two model self-assembled QWR structures, 
one with (2/2) SPS structure and the other with (2/2.25) SPS structure, and with varying degrees 
of lateral alloy mixing. The clamping effect due to thick capping layer is also examined. The valence 
force field (VFF) model is used to calculate the equilibrium atomic positions in the model QWR 
structures. This allows the calculation of the strain distribution at the atomistic level. 

We found that in model structures with thick capping layer (clamped case), the magnitude of 
shear strain is quite large (around 2%) for abrupt lateral composition modulation and it reduces 
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substantially when the lateral alloy mixing is introduced. We found that the shear strain effect can 
alter the valence subband structures substantially when the hole is confined in the Ga-rich region, 
and it gives rise to a reversed optical anisotropy compared with QWR structures with negligible shear 
strain. This points to a possibility of " shear strain-engineering" to obtain QWR laser structures of 
desired optical anisotropy. For model structures with thin capping layer (undamped case), the effect 
of shear strain is less dramatic, since the hole is confined in the In-rich region, where both the QWR 
confinement and the shear strain have similar effect. We found that polarization ratio (-Pfii ]/-P[iio]) 
for the undamped case is always larger than 1 for both (2/2) and (2/2.25) SPS structures, regardless 
of the degree of lateral alloy mixing. 

The band gap and polarization ratio obtained for both clamped and undamped (2/2) SPS struc- 
tures with lateral alloy mixing (x rn ps 0.8) and for undamped (2/2.25) SPS structures are consistent 
with experimental observations on most self-assembled InGaAs QWRs with similar specifications. 
For clamped (2/2) SPS structure without alloy mixing, we predicted a reversed optical anisotropy 
compared with the PL measurements for most InAs/GaAs self-assembled QWRs. However, there 
exist a few InAs/GaAs self-assembled QWR samples which display the reversed optical anisotropy 
at 77 K[7], and recent experimental studies also indicated a reversed optical anisotropy in many 
samples at 300 K.[6] Our studies demonstrated one possibility for the reversed optical anisotropy, 
the shear strain effect. There are other possibilities for the reversed optical anisotropy at 300K, such 
as the thermal population of the excited hole states, which can have different atomic characters as 
discussed in the previous section. The temperature dependence of the PL peak and polarization in 
self-assembled InGaAs QWRs[6] remains an intriguing question to be addressed in the future. 
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Table I. 



List of band gaps and polarization ratios for all structures studied 



x m b/ano alloy mixing clamped/unclamped SPS structure gap(eV) P\\/P± 



1.0 





no 


clamped 


(2/2) 


0.765 


0.36 


0.8 


7 


yes 


clamped 


(2/2) 


0.763 


3.5 


1.0 





no 


undamped 


(2/2) 


0.791 


3.1 


0.7 


7 


yes 


undamped 


(2/2) 


0.766 


1.8 


0.6 


7 


yes 


undamped 


(2/2) 


0.80 


1.82 


0.6 


15 


yes 


undamped 


(2/2) 


0.798 


1.5 


1.0 





no 


undamped 


(2/2.25) 


0.731 


2.2 


1.0 


7 


yes 


undamped 


(2/2.25) 


0.74 


2.1 
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Figure Captions 

Fig. 1. Schematic sketch of the unit cell of the SILO quantum wire for two model structures con- 
sidered. Each unit cell consists of 8 pairs of (2/2) or (2/2.25) GaAs/InAs short-period superlattices 
(SPS). In structure 1, the (2/2) SPS contains four diatomic layers (only the cation planes are il- 
lustrated) and the SPS is repeated eight times. In structure 2, four pairs of (2/2.25) SPS (or 17 
diatomic layers) form a period, and the period is repeated twice in the unit cell. Filled and open 
circles indicate Ga and In rows (each row extends infinitely along the [110] direction). 
Fig. 2. Diagonal strain distribution of model structure 1 [corresponding to Fig. 1(a)] with thick 
capping layer (clamped case). Solid: e x / x /. Dashed: e y i y i. Dash-dotted: e zz . Parameters for alloy 
mixing: x m = 0.8 and b = 7a [no]. 

Fig. 3. Diagonal strain distribution of model structure 1 with thin capping layer (undamped case). 
Solid: e x > x i. Dashed: e y > y >. Dash-dotted: e z > z >. Parameters for alloy mixing: x m = 0.8 and b = 7a[no]. 
Fig. 4. Diagonal strain distribution of model structure 2 [corresponding to Fig. 1(b)] with thin 
capping layer (undamped case). Solid: e x i x t. Dashed: e y i y >. Dash-dotted: e z > z i. Parameters for alloy 
mixing: x rn = 0.8 and b = 7a[n ]. 

Fig. 5. Conduction and valence band edges for 8 pairs of (2/2) SPS structures sandwiched between 
6OAAlo.24Gao.24Ino.52As barriers with different degrees of lateral alloy mixing as functions of the [110] 
coordinate y' for both clamped (thick capping layer) and undamped (thin capping layer) case. In 
upper pane, the dashed line is for abrupt modulation (without alloy mixing) and the solid line is for 
lateral alloy mixing with x m = 0.8 and b = 7a [no]. In lower panel, dash-dotted line: without alloy 
mixing; solid line: x m = 0.7 and b = 7a[noj; dotted lines: x m = 0.6 and b = 7a[no]; dashed lines: 
x m = 0.6 and b = 15ariio]. 

Fig. 6. Conduction and valence band edges for 8 pairs of (2/2.25) SPS structures sandwiched 
between 6OAAlo.24Gao.24Ino.52As barriers with different degrees of lateral alloy mixing as functions of 
the [110] coordinate y' for both clamped (thick capping layer) and undamped (thin capping layer) 
case. Dashed lines: without alloy mixing. Solid lines: with alloy mixing described by x m — 1.0 and 

b = 7a[no]. 

Fig. 7. Valance subband structures and squared optical matrix elements for self-assembled QWR 
made of (2/2) SPS structure with thick capping layer (clamped case) and without alloy mixing. 
Fig. 8. Valance subband structures and squared optical matrix elements for self-assembled QWR 
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made of (2/2) SPS structure with thin capping layer (undamped case) and with alloy mixing described 
by x m = 0.8 and b = 7a [no]. 
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Table. 1 



The list of results of all structures studied 



x m b(ano) alloy mixing clamped/unclamped structure model gap(eV) P\\/P±_ 

1.0 no clamped 1 0.765 0.36 

0.8 7 yes clamped 1 0.763 3.5 

1.0 no undamped 1 0.791 3.1 

0.7 7 yes undamped 1 0.766 1.8 

0.6 7 yes undamped 1 0.80 1.82 

0.6 15 yes undamped 1 0.798 1.5 

1.0 no undamped 2 0.731 2.2 

1.0 7 yes undamped 2 0.74 2.1 
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